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ABSTRACT 

We present positions, kinematics, and the planetary nebula luminosity function 
(PNLF) for 35 planetary nebulae (PNe) in the nearest starburst galaxy ICIO extend- 
ing out to 3 kpc from the galaxy's centre. We take advantage of the deep imaging and 
spectroscopic capabilities provided by the spectrograph FOCAS on the 8.2m Subaru 
telescope. The FN velocities were measured through the slitless-spectroscopy tech- 
nique, which allows us to explore the kinematics of ICIO with high precision. Using 
these velocities, we conclude that there is a kinematic connection between the H i 
envelope located around ICIO and the galaxy's FN population. By assuming that the 
FNe in the central regions and in the outskirts have similar ages, our results put 
strong observational constraints on the past tidal interactions in the Local Group. 
This is so because by dating the FN central stars, we, therefore, infer the epoch of a 
major episode of star formation likely linked to the first encounter of the H i extended 
envelope with the galaxy. Our deep [O ill] images also allow us to use the FNLF to 
estimate a distance modulus of 24.1±0.25, which is in agreement with recent results 
in the literature based on other techniques. 

Key words: galaxies: Local Group; galaxies: individual: ICIO; galaxies: kinematics 
and dynamics; ISM: planetary nebulae 



1 INTRODUCTION 

Even though dwarf galaxies are currently the most abundant 
galaxy population, and were probably even more abundant 



formation is still under debate (|Matedll998 


). According to 


the hierarchical formation scenario (see, e.g.. 


White & FrenkI 



so their study is crucial for understanding galaxy formation 
and evolution processes. One of the several questions about 
their formation and evolution, and that still remains un- 
solved, is the origin of the differences be tween the irr egular 
and spheroidal morphological types fsee lMateolll998i . for a 
review) . We have been approaching this issue by using plane- 
tary nebulae (PNe, old-intermediate stellar population) and 
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H II regions (young population) as tracers of the chemi- 
cal evolution of the Local Group (LG) dwarf galaxies (see, 
for example, Goncalves et al.l 20071: Magrini fc Goncalved 
l2009l : iMagrini et al.l I2OIII : iGoncalves et al.l \20li Y In par- 
ticular, we have dedicated the last few years in the study 
of the d warf irregular galaxy ICIO - its PNe and H 11 
regions (|Magrini et all l2003l : iMagrini fc Goncalved 120091 '): 
its serendipitously discovered symbiotic s ystem, which h 
the farthest known symbiotic star to date (Goncalve s et al 
|2008); and its detailed chemical evolution history (,Yin et al 
l2010f ). Here we present a kinematic study of PNe in this 
galaxy. As a by-product, we also present some photometric 
results. 

ICIO is a dwarf irregular galaxy located in the outskirts 
of the Local Group at a low galactic latitude (1 = 119°. 0, b 
= -3°. 3; NASA/IPAC Extragalactic Database, NED), which 
makes the reddening values very uncertain. The estimates 
of the reddening range from E(B-V)=0.47 to 2.0, translat- 
ing into a distance modulus (m — Al) that varies from 22 to 
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more than 27 mag, i.e., the linear distance for IC 10 would 
vary from 0.5 to 3 Mpc (Sakai, Madore & Freedman 1999; 
Demers, Battinelli & Letarte 2004). More recent estimates 
locate ICIO at a distance around 0.7-0.8 Mpc (e.g. , Pustilnik 
& Zucker 2008: ISanna et al.ll2008l : iKim et alll2009l ). suggest- 
ing a possible membership to the M31 subgroup. 

In addition, ICIO has a high star formation rate 

i SFR) as evidenced by its large number of H 1 1 regions 
Hodge fc Le3 Il990l) as well as fo r its Ha (jMateol [l998l ) 
and far-IR ( Melisse fc Israeli 1 1994 ) luminosity. It also has 
a huge number of Wolf-Rayet stars per unit luminosity, 
the largest in the LG (Massey, Armandroff & Conti 1992; 
iMassev fc Armandrofj 1 19951 ). and an anomalous ratio of 
carbon-type Wolf-Rayet stars (WC stars) to nitrogen-type 
Wolf-Rayet stars (WN s tars), which is pecul iar at its metal- 
licit y (0.2-0.3 sola r; cf. ISkillman. Kennicutt fc Hodgdll989l 
and ' lGarnett|[l990l)^ The presence of a.n ext ended and com- 
plex H I envelope jshostak fc Skillmanll 19891 . hereafter SS89) 
suggests that ICIO is still in its formative stage or that it has 
undergone a major merger episode ( Wilc ots fc Mil ler 1998, 
hereafter WM98). Altogether these characteristics suggest 
that ICIO is experiencing an intense and very recent burst 
of star formation, likely triggered by infalling gas from its 
H I envelope which is counter-rotating with respect to the 

galaxy proper rotation (WM98). 

The Local Group Census (cf. ICorradi fc Magrinill2006l ) 
discovered a large number of planetary nebulae (PNe) in 
many dwarf galaxies, including ICIO. Several of these PNe 
were found at large distances from the centre of their host 
galaxies, and so could represent a trace of past tidal interac- 
tions and/or a connection of the old stellar population with 
the H I envelope located around some dwarf galaxies, as in 
the case of ICIO (Huchtmcior 1979). 

Single-dish 21 cm observations showed that ICIO is em- 
bedded within a huge envelope of H I extending seven times 
farther than its o ptical diameter an d stretching nearly a de- 
gree on the sky l|Huchtmeieij 1 197*91 : Huchtmeier, Seiradakis 
fc Materne 1981). Moreover, the H I velocity field of ICIO 
appears to be quite complex (|Huchtmeierlll979l ) with a long 
stream extending to the South of the galaxy. Interferometric 
data revealed that ICIO has a rotating disk surrounded by a 
counter-rotating or highly warped outer distribution of gas 
(SS89). 

The most likely scenario is that the gas component rep- 
resents a late acquisition. Numerical simulations, for exam- 
ple, i ndicate that merger s can lead to counter-rotating cores 
(e.g., iThakar et al]|l997l ). The detection of a stellar popu- 
lation related to the outer envelope would be extremely im- 
portant to constrain its age and thus the evolution of ICIO, 
dating the possible encounter between the two components. 
However, so far deep optical images did not show any stel- 
lar population coincident with the extended gas envelope, 
as discussed by WM98 from their examination of Digital 
Sky Survey and deep WIYN R band images. The surface 
brightness of this stellar population is expected to be very 
low, and its detection is strongly limited by scattered light 
from nearby bright objects and by the contribution of faint 
discrete sources. These limitations can be circumvented ob- 
serving PNe because they are easily detectable from their 
strong nebular emission of [O ill] at 5007 A. 

In this paper we aim at studying the radial velocities of 
the PN population of ICIO in order to explore the connec- 



tion of its intermediate to old stellar population with its H I 
envelope. Previous to the present work, the PN kinemat- 
ics of this galaxy was never explored, although a number 
of PNe were already known. Altogether sixteen candidate 
P Ne were identi fi ed in ICIO by the Local Group Census 
m iMagrini et all (|2003l ). Other nine fainter PNe were de- 
te cted and immediately c onfirm ed using deep spectroscopy 
bv lMagrini fc Goncalvej (|2009| ). 



2 DATA ACQUISITION AND REDUCTION 

The observations were made by one of us (R.H.M.) with 
the Faint Object Ca mera and Spectrograph (FOCAS; 
iKashikawa et al.ll200d ) attached to the Cassegrain focus of 
the 8.2 m Subaru telescope, Mauna Kea, Hawaii. 



2.1 Observations 

Data were obtained during a single night, 2010 September 
30th. The night was photometric, with an average seeing of 
0".8. 

The field of view of FOCAS is 6.5 arcmin and is cov- 
ered by two CCDs of 2k x 4k (pixel size 15 ^m) with an 
unexposed gap of 5" between them. For simplicity, we will 
call the two CCDs Chip 1 and Chip 2. The image scale is 
0.104 arcsec pixel"^. 

The observations were obtained with two filters: the on- 
band filter with a central wavelength of 5025 A, a FWHM 
(Full Width Half Maximum) of 60 A, a peak transmission 
of 68%, and an equivalent width of 40 A, and the off-band 
filter image was done through a standard broadband visual 
filter. We observed three fields - A, B, and C - centred re- 
spectively at {a,S, J2000) = (Oh 19m 00.4s, 59°16'25"; Oh 
20m 03.9s, 59°17'18"; Oh 20m 56.9s, 59°16'48"), as shown 
in Figure [1] The observing sequence for each field was: one 
off-band image (exposure time 150 s), one on-band image 
(exposure time 1200 s) and two grism-|-on-band images (ex- 
posure time 1200 s). This observing sequence was executed 
twice for fields A and B and three times for field C. 

For the wavelength calibration, after each grism-|-on- 
band exposure, an engineering mask was inserted in the 
light path, and on-band and grism-f on-band images were 
obtained illuminating the mask with the comparison lamp. 
In addition, on-band and grism-|-on-band images of the en- 
gineering mask were taken illuminating the mask with the 
Galactic PN NGC 7293 (PNG 036.1-57.1), for radial veloc- 
ity quality control. E xamples of the calibration images are 
shown in Figure 2 of iMendez et al.l l|2009l) . Table 1 shows 
the log of the observations with the corresponding exposure 
times and air masses. For brevity, we list only the images 
corresponding to Chip 1. Identification numbers of Chip 2 
images are obtained by adding 1 to the Chip 1 identification 
number s in Table 1. The spectrophotometric standard G 
138-31 (|Okell 19901 ) was used for the photometric calibration 
of the on-band images. 



2.2 (Re)identiflcation of PNe, FOCAS slitless 
spectroscopy and Reductions 

The traditional on-band/off-band filter technique was used 
for the re-detection of the known PNe and for the iden- 
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Figure 1. The distribution of all PN candidates (squares) de- 
tected with Subaru+FOCAS in ICIO. The big circles indicate the 
6.5 arcmin FOCAS field of view. At the distance of 660 kpc (see 
section 4.2), one arcmin corresponds to 0.19 kpc. 



Table 1. Journal of observations for Chip 1. 



FOCAS Field Configuration FOCAS id Exp Air Mass 



G 138 


- 31 




On-band 


120371 


60 


1 


.33 


G 138 


- 31 




On-band 


120373 


60 


1 


.36 


NGC 7293 + mask 


On-band 


120397 


300 


1 


.87 


NGC 7293 + mask 


On-t-grism 


120399 


600 


1 


.80 


IC 10 


Field 


B 


Off-band 


120403 


150 


1 


.94 


IC 10 


Field 


B 


On-band 


120405 


1200 


1 


.86 


IC 10 


Field 


B 


On-t-grism 


120407 


1200 


1 


.75 


IC 10 


Field 


B 


On-t-grism 


120409 


1200 


1 


.65 


IC 10 


Field 


B 


Off-band 


120423 


150 


1 


.46 


IC 10 


Field 


B 


On-band 


120425 


1200 


1 


.43 


IC 10 


Field 


B 


On-t-grism 


120427 


1200 


1 


.39 


IC 10 


Field 


B 


On-t-grism 


120429 


1200 


1 


.36 


IC 10 


Field 


A 


Off-band 


120431 


150 


1 


.33 


IC 10 


Field 


A 


On-band 


120433 


1200 


1 


.32 


IC 10 


Field 


A 


On-t-grism 


120435 


1200 


1 


.31 


IC 10 


Field 


A 


On-t-grism 


120437 


1200 


1 


.30 


IC 10 


Field 


C 


Off-band 


120443 


150 


1 


.30 


IC 10 


Field 


C 


On-band 


120445 


1200 


1 


.30 


IC 10 


Field 


C 


On-t-grism 


120447 


1200 


1 


.31 


IC 10 


Field 


C 


On-t-grism 


120449 


1200 


1 


.32 


IC 10 


Field 


C 


Off-band 


120451 


150 


1 


.33 


IC 10 


Field 


C 


On-band 


120453 


1200 


1 


.35 


IC 10 


Field 


C 


On-t-grism 


120455 


1200 


1 


.38 


IC 10 


Field 


C 


On-t-grism 


120457 


1200 


1 


.42 


IC 10 


Field 


A 


Off-band 


120459 


150 


1 


.45 


IC 10 


Field 


A 


On-band 


120461 


1200 


1 


.48 


IC 10 


Field 


A 


On-t-grism 


120463 


1200 


1 


.54 


IC 10 


Field 


A 


On-t-grism 


120465 


1200 


1 


.62 


IC 10 


Field 


C 


Off-band 


120471 


150 


1 


.70 


IC 10 


Field 


C 


On-band 


120473 


1200 


1 


.76 


IC 10 


Field 


C 


On-t-grism 


120475 


1200 


1 


.88 


IC 10 


Field 


C 


On-t-grism 


120477 


1200 


2 


.03 



tification of new PN candidates. The on-band image is 
taken through a narrow-band filter passing the redshifted [O 
III] 5007 A nebular emission line, while the ofT-band image is 
taken through a broader filter passing no nebular emissions. 
The PNe are visible as point sources in the on-band image, 
but are absent in the ofT-band image. A third image, taken 
through the on-band filter and a grism, confirms the PN can- 



Table 2. Radial velocities and magnitudes for the detected 
PNe in ICIO. 



ID a (5 RV mpoCAS mAf03 

(.12000) {J2000) kms-1 (5007 A) (5007 A) 



01-A 





18 


44 


,61 


59 


17 


47 


,83 


-328 


23 


,83 


24.12 (1) 


02-A 





19 


04 


,40 


59 


17 


04 


,82 


-350 


23 


,63 


24.01 (2) 


03-A 





19 


13 


,64 


59 


14 


44 


,55 


-361 


23 


,14 


23.58 (3) 


04-B 





19 


49 


,37 


59 


17 


11 


,08 


-371 


24 


,83 




05-B 





19 


57, 


,71 


59 


20 


31 


,42 


" 


26 


,52 




06-B* 





19 


57 


,51 


59 


17 


11 


,37 


-317 


24 


,71 




07-B 





19 


59 


,37 


59 


18 


42 


,11 


-361 


25 


,00 




08-B 





20 


01 


,15 


59 


16 


41 


.55 


-365 


24 


,87 




09-B 





20 


03 


,65 


59 


20 


44 


,05 


" 


24 


,71 




10-B* 





20 


03 


,90 


59 


19 


26, 


,12 


-323 


24 


.65 




11-B 





20 


04 


,20 


59 


16 


55 


,30 


-340 


24 


.18 




12-B 





20 


06 


,34 


59 


17 


05 


,24 


-308 


25 


,89 




13-B 





20 


06 


,66 


59 


15 


22 


,18 


-354 


23 


,94 


24.76 (4) 


14-B 





20 


09 


,68 


59 


19 


47 


,17 


-310 


23 


,83 




15-B 





20 


09, 


,87 


59 


16 


34, 


,67 


-357 


25 


,10 




16-B 





20 


10, 


,79 


59 


19 


56, 


,67 


-344 


23 


,92 


- 


17-B 





20 


12 


,75 


59 


19 


11 


.53 


-313 


25 


,73 




18-B* 





20 


14 


,74 


59 


18 


08 


.14 


-366 


25 


,60 




19-B 





20 


14 


,89 


59 


18 


53 


,46 


-320 


24 


,50 




20-B 





20 


15 


,82 


59 


14 


59 


,75 


-333 


26 


,19 




21-B 





20 


16 


,21 


59 


17 


07 


,33 


-377 


26 


,22 




22-B 





20 


17 


,03 


59 


17 


58 


,24 


-327 


24 


,67 




23-B* 





20 


17 


,17 


59 


15 


52 


,71 


-342 


23 


,64 


24.49 (5) 


24-B 





20 


17, 


,87 


59 


18 


08, 


,72 


-365 


25, 


,28 




25-B 





20 


18, 


,00 


59 


18 


11 


,48 


-333 


25, 


,69 




26-B 





20 


18, 


,14 


59 


18 


58, 


.17 


-285 


25, 


,27 




27-B 





20 


19, 


,32 


59 


18 


02, 


,70 


-324 


22 


,54 




28-B* 





20 


19, 


,66 


59 


18 


13, 


,03 


-341 


24, 


,02 




29-B 





20 


22 


,28 


59 


16 


46, 


,59 


-340 


25, 


,23 




30-B 





20 


24, 


,28 


59 


15 


38, 


,81 




25, 


,41 




31-B 





20 


27, 


,58 


59 


16 


53, 


,83 




26, 


,12 




32-B 





20 


27, 


,98 


59 


18 


20, 


,91 




24, 


,78 




33-C 





20 


40, 


,65 


59 


17 


46, 


,32 


-284 








34-C 





20 


40, 


,89 


59 


18 


50, 


,15 




22 


,93 


23.33 (10) 


35-C 





20 


42 


,95 


59 


16 


30, 


,90 


-332 


22 


,85 


23.43 (11) 


36-C 





20 


46, 


,34 


59 


18 


38, 


,30 


-372 


26, 


,43 




37-C 





20 


47, 


,40 


59 


19 


00, 


,05 


-284 








38-C 





20 


49, 


,06 


59 


18 


33, 


,16 


-279 


24, 


,18 




39-C 





21 


09, 


,80 


59 


16 


59, 


,70 


-331 


22 


,81 


23.18 (15) 


40-C 





21 


11 


,20 


59 


16 


11 


,71 


-382 


23, 


,81 


24.52 (16) 


41-C 





21 


14, 


,40 


59 


18 


08, 


17 


-303 


26, 


,37 





Notes — (i) The asterisks in the ID column indicate the confirmed 
PNe studied, through deep medium-resolution optical spectroscopy, 
by iMagrini &: Gongalved j2009|) . (ii) The last two columns present the 
magnitudes obtained in this work (mpOCAS! refer to the text for more 
details) and in lMagrini ct al.l |2003l , M03, in which the photometry and 
astrometry of a number of PN candidates in ICIO were discussed. See 
their Tabic 1). In the last column, the numbers in parentheses stand 
for the M03's IDs. 

didates. By inserting the grism in the light path, the images 
of all continuum sources are transformed into segments of 
width determined by the on-band filter transmission curve. 
All the emission-line point sources, such as the PNe, remain 
as point sources. The grism also introduces a shift relative 
to the undispersed on-band image which is a function of the 
wavelength of the nebular emission line and of position on 
the CCD. By calibrating this shift, we are able to measure 
the radial velocities for all emission-line objects in the field, 
as in Mendez et al. (2001, 2009). In the case of Subaru and 
FOCAS, the dispersing element was an echelle grism with 
175 grooves/mm which operates in the 4th order and gives 
a dispersion of 0.5 A pixel"^, with an efficiency of 60% (see 
the Subaru website). 

Standard irafQ tasks were used for the basic CCD re- 
ductions (bias subtraction, flat field correction using twilight 

^ IRAF is distributed by the National Optical Astronomical Ob- 
servatories, operated by the Association of Universities for Re- 
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Average Velocity 
335±28 km s-i 



-350 -300 
Heliocentric RV (km s"') 



Figure 2. Histogram of heliocentric radial velocities for 35 PN 
candidates in ICIO, in bins of 20 kms~^. The average velocity 
and dispersion are given in the figure. 



flats). We refer the reader to the papers by Mendez et al. 
(2001, 2009) where a detailed description of all the steps 
involved in the reduction analysis is given. 



3 KINEMATIC ANALYSIS 

In Figure [1] we show the spatial distribution of the 
41 PN candidates that we have identified in ICIO with 
Subaru+FOCAS. As explained in Sect. 2.2, we have used 
the slitless technique to obtain the heliocentric radial ve- 
locities (RVs) for the sample. Hereafter, we will refer to 
our RVs simply as 'velocities'. Typical images with seeing 
around 0.5" have a point-spread function (PSF) size of 5 
pixels, which translates into a radial velocity resolution of 
140 km s~^, i.e., the internal velocity field of the PN is not 
resolved. Assuming position errors of 0.4 pixel, the expected 
uncertainty in radial velocity is 12 km s~^. However, due 
to background contamination, bad columns in the images, 
and the displacement produced by the grism placing some 
of the sources outside of the fields, we were able to measure 
the velocities of 35 objects only. The final values are listed in 
Table 2. Adding quadratically the uncertainties in the veloc- 
ities, which are given by the calibration errors (~10 kms~^), 
the position errors (~10 kms~^), and the errors from image 
registration (~10 kms~^), we get an overall error of about 
17 kms~^. Assuming that the spectrograph deformations 
and guiding errors have a marginal contribution, we esti- 
mate the total uncertainty in the velocities measured with 
Subaru to be at most 20 kms~^. 

We show in Figure [2] a histogram of the velocities mea- 
sured for our targets. There is a clear velocity peak at 
— 335kms~^, with a dispersion of ±28 kms"\ Within the 
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uncertainties, the average PN velocity is in good agreement 
with the NED heliocentric radial velocity of — 348kms~^ 
for ICIO. We note from the velocity histogram that none 
of the objects is apart from the others in terms of velocity 
distribution. This finding suggests that the objects are true 
members of ICIO. 

We then obtained the galactocentric coordinates x and 
y for the objects. For that purpose we have used a central 
position of a(J2000)=0h 20m 17.3s and (5(J2000)=59° 18'14" 
for ICIO. As ICIO is an irregular galaxy, i t is dif ficult to de- 
fine its dynamical axis. IWilcots fc Milleij { 19981 ) argue that 
even though somewhat uncertain the position angle (PA) 
value might range between 50 and 90 degrees (see. e.g., their 
Figure 5). For this reason, we chose to investigate any sign 
of rotation coming from the PNe along axes with different 
PA, slightly extending the WM98's PA range. Along each 
of those axes we defined four subsamples and calculated the 
average velocity for each one. We note, however, that the 
subsamples were not equally spaced in x because we want 
to ensure a minimum number of PNe per subsample. Fig- 
ures 3 to 6 show the binned velocities of the PNe along axes 
with PA of 50°, 60°, 90° and 135°, respectively. We selected 
these angles because PA= 135° corresponds to the major axis 
of the IC 10 starlight distribution (as can be seen in Figure 
1), while PA=60°gives the direction along which the neutral 
gas is rotating, according to SS89 and WM98. 

Figures [3] |4] and [5] suggest some rotation along the axis 
with PA= 50°, 60°and 90°, while figure [6] shows that the 
evidence of rotation along the optical PA= 135° axis is less 
significant. We conclude that the PNe (and therefore the 
stellar population) appear to be rotating like the neutral 
gas. 

Let us make a more detailed comparison. In WM98, 
the authors produced velocity maps at different locations 
in the sky, each map corresponding to a velocity channel. 
In Figure [7] the positio n of our PNe wit h velo cities falling 
within the range of the IWilcots fc Milled (|l998!) map chan- 
nels are superimposed to the corresponding map. We note 
that in most maps PNe are found to coincide with the gas, 
further suggesting that the PN kinematics follows the H I 
one. Note that WM98 concluded, using the velocity maps we 
just discussed, that their H I velocities followed the rotation 
curve obtained previously by SS89. So we can also compare 
our PN velocities to the plots in SS89. In Figure 8 we have 
taken the position-velocity diagram from SS89 (their Fig. 
10) and we have superimposed our PN velocities. This dia- 
gram was constructed with an x coordinate measured along 
PA=60°, following the observed neutral gas rotation. Along 
this direction, the H I velocities are very well behaved, as is 
most clearly seen in SS89's Fig 8. Turning to our objects, in 
spite of the high PN velocity dispersion, it is immediately 
apparent that the PNe are rotating in the same direction, 
with less negative velocities in the upper left. There is a 
good overlap with the neutral gas velocities, although sev- 
eral PN velocities are less negative. The full interpretation 
of this figure is not straightforward, because, although the 
X coordinates are the same, the y coordinates can be quite 
different. Consider for example the three PNe in Field A 
(see Figure 1). Their x coordinates are the most negative 
in our sample (< —500), and they might be expected to lie 
at the most negative velocities (they appear at the extreme 
right of Figure 8); but their average velocity is —340 km s~^. 
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Figure 3. Left: the same spatial distribution of PNe as in Figure [T] The dashed Hne represents an axis of PA = 50°. To study the 
rotation along this axis, the PN sample was divided into four bins perpendicular to it. The orientation of this axis is shown on the top 
and the PNe inside each bin are represented as: open circles, filled triangles, open squares and stars. Right: the average velocity within 
each bin is plotted at the position of the average x-coordinate of the bin. The 4 bin ranges are: x < —100" (open circle), —100" < a; < 0" 
(filled triangle), 0" < x < 100" (open square), x > 100" (stars). 
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Figure 4. As for Figure [S] except that the rotation is studied along an axis of PA=60°. 
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Figure 5. As for Figure [S] except that the rotation is studied along an axis of PA=90° 
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Figure 6. As for FigureO except that the rotation is studied along an axis of PA=135° 
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Figure 8. The position-velocity diagram from Figure 10 of SS89, 
showing the kinematics of neutral gas along PA=60°. The squares 
indicate our PN velocities as a function of the corresponding x 
coordinate, calculated along an axis with PA=60°, as in Figure 
4. 



4 PN PHOTOMETRY 

A detailed d escription of the PN pho tometric procedure is 
provided by iTeodorescu et aEI (|201Cl ). The results of our 
photometry for 39 of the 41 objects are listed in Table 2. 
Two objects (PN 33 and 37) were too close to bright fore- 
ground stars, which precluded accurate photometry. The [O 
III] 5007 A fluxes measured through the on-band filter are 
traditionally expressed in magnitudes m(5007), using the 
definition introduced by Jacoby (1989,). 

m(5007) = -2.5 log /(5007) - 13.74. (1) 

For the flux ca libration, we adopted the standard star G138- 
31 (|Okelll99(]| '). which has a monochromatic flux at 5028 A 
(the redshifted wavelength of its [O ill] 5007A emission due 
to the PN radial velocity), of 1.44 x 10"^^ ergs cm~^ s~^ 

A-^ 

If we compa re our photometry with that by 
iMagrini et al.l (|2003l ) for the objects in common, we flnd that 
our values are systematically brighter, by abo ut 0.4 mag. We 
attribute this to the fact that the nights of iMagrini et al.l 
(2003) were not of photometric quality, and thus our new 
values presented in Table 2 are preferred. 



However, they should not be compared with the neutral gas 
velocity distribution shown in Figure 8, because they are in 
a different part of the sky, to the northwest. In fact, there 
is a blob of H i close to these Field A PNe, seen in SS89's 
Fig. 9, and it has a velocity of about —330 km s~^, which 
flts those few PN velocities quite well, perhaps by chance. 

In summary, leaving aside some discrepancies, which 
probably but not necessarily indicate a tendency toward 
more negative neutral gas velocity than PN velocities, in 
Figs. 3 to 8 we have been able to show that most of the PN 
velocities follow closely the H l kinematics pattern derived 
by SS89 and WM98. This illustrates the regular, solid-body 
rotation-like nature of the velocity fleld (given by the H I as 
weU the PNe) in IClO's disk. 



4.1 Foreground extinction 

Being located at a low galactic latitude (1=119.0°, b=-3.3°), 
ICIO is affected by signiflcant foreground extinction. As 
previously stated, we flnd in the literature different v alues 
for this parameter. ISchlegel. Finkbeiner fc David (| 19981 ) flnd 
E[B — V) = 1.527, while studies based either o n the tip of 
the r ed giant branch ( TRGB:,Sakai et al. 1999. : .Sanna et al.l 
I2008D or the Cepheids (| Wilson et al.ll 19961 7 give a range be- 
tween 0. 6 to 1.1 mag. We h ave adopted E{B -V) = 0.77, 
following [Richer et al.l (|200ll ). whose extinction value is pre- 
ferred because of the very convincing work they did, by plot- 
ting the reddening as a function of the H I column density 
in their Figure 3. 
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Figure 7. The PN (filled circles) positions superimposed on the H i maps produced by WM98. The velocities of the PNe in each map 
fall within the velocity range of the map. 
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Figure 9. Observed [O III] 5007 PNLF of ICIO (squares), with the 
sample of 35 PNe binned into 1-mag intervals. The apparent mag- 
nitudes m(5007) have been transformed into absolute magnitudes 
M(5007) by adopting an extinction correction of 2.74 mag and a 
distance modulus of 24.1. The three lines are PNLF simulations 
(Mendez and Soffner 1997) for 3 different total PN population 
sizes: 100, 150 and 200. 



4.2 The distance determination from the PN 
luminosity function 

Having measured the apparent magnitudes m(5007) of the 
35 PNe, the PN luminosity function (PNLF) was built, using 
1-mag b ins, and compared with s imulated PNLFs like those 
used bv lMendez fc Soffneij |l993) to fit the observed PNLF 
of M31. The comparison is shown in Figure |9l The absolute 
magnitudes M(5007) that produce the best fit to the simu- 
lated PNLF were calculated usin g an extinction corre ction 
A(A), at 5007 A, of 3.56£'(B - V) (jCardelh et al.lll989l ). and 
a distance modulus (m — M) = 24.1, which is equivalent to 
660 kpc. 

The simulated PNLFs plotted in Figure [9] are binned, 
like the observed one, into 1-mag intervals and have a max- 
imum final mass of 0.63 Ma and sample sizes between 100 
and 200 PNe fsee lMendez fc SoffneiHToOTl : the "sample size" 
is the total number of PNe, detected or not, that exist in 
the surveyed area). We selected a rather large bin size of 
1 mag to compensate for the rather small number of PNe, 
and reduce the statistical noise in the PNLF. We estimated 
an error of 0.2 mag from the goodness of the fit at differ- 
ent distance moduli. To obtain the total error estimate, we 
have to combine the possible systemat ic and random errors . 
The systematic error is the same as in Ijacobv et aP (Il990l i. 
i.e., 0.13 mag, including the possible error in the distance to 
M31, in the modeling of the PNLF and in the foreground 
extinction. The random contributions are given by 0.2 mag 
from the fit to the PNLF, 0.05 mag from the photometric 
zero point, and 0.05 mag from the filter calibration. Combin- 
ing all these errors quadratically, we estimate that the total 
error bar for the PNLF d istance modulu s is ± 0.25 mag. 
A PNLF was obtain ed by iKniazev et al.l (|2008l ') , based on 
iMagrini et al.l l|2003h magnitudes. Given th e uncertainties in 
both determinations. iKniazev et al.l (|2008l Vs distance mod- 



ulus (24.30_(] jg) is in good agreement with o urs. Note that 
we ha ve adopted bri ghter magnitu d es th an iMagrini et al.l 
(':2003l '). If we correct IKniazev et al.l (|2008l Vs distance mod- 
ulus by the same amount, their distance modulus decreases 
to 23.9. 

By studying the TRGB in ICIO ISakai et al.l (|l999l ') 
found a distance modulus of 24.1±0.2 mag fo r Population 
I and 23.5±0.2 mag for Population II, while I S anna et al.l 
(: 2008l l estimated a distance modulus of 24.57±0.15 mag. 
From the observations of the Cepheids. I Wilson et all (| 19961 ) 
deduced a distance modulus of 24.57±0.21 mag. Note that 
there may be some amount of internal extinction in ICIO, 
in addition to the foreground extinction. For the sake of a 
complete knowledge of the ICIO properties, although not 
significantly affecting our conclusions about the kinematic 
pattern of the galaxy, it would be important to know how 
much extra extinction to expect. Obviously we cannot go 
further on this discussion without deep individual spectra 
for all the PNe we found. Such a work, clearly, lies outside 
of the scope of the present paper. Had we increased the ex- 
tinction correction t o 3.5 mag, equivalent to E{B — V) = 1.0 
(jPemers et al.l [20041 '). our new distance modulus would be- 
come 23.3. Such a large correction is unlikely, because many 
of the objects we discovered are in fact in the outskirts of 
ICIO. We give these numbers to emphasize that the value of 
the extinction correction is by far the most important source 
of uncertainty in the PNLF distance determination for ICIO. 

4.3 PN formation rate 

Once the sample size is known, we can calculate the specific 
PN formation rate ^ in units of PNe yr~^ Lq~^: 

npN — (.LrtpN, (2) 

where npN is the sample size, Lt is the total bolometric 
luminosity of the sampled population, expressed in Lq , and 
tpN is the lifetime of a PN, for which we have adopted 
30,000 yr in the PNLF simulations. We have the B ban d 
mag, Bt=8.68, and B - V=0.53, from [Richer et al.l (1200 jj. 
and a bolometric correction of —0.8 mag ( Buzzoni et al.l 
1 2006) from which we obtain an extinction corrected appar- 
ent bolometric magnitude 7.35. Using the distance modulus 
we obtained, of 24.1, and a solar Mboi = 4.72 we calculate 
the total luminosity of IC 10, Lt = 3.9 x 10** Lq. The sam- 
pled luminosity is 90% of the total luminosity, i.e. 3.5 x 10* 
Lq. Adopting npM = 150 (see subsection above), we get ^ = 
(14.28 ± 5) X 10"^^ Transforming our value of ^ into a spe- 
cific PN density -or a ratio, defined as the number of PNe 
to the total galaxy luminosity ratio-, as log (a) = —6.37, in 
agreement with IClO'starburst quality (|Buzzoni et al.ll2006l 
obtained —6.59). 



5 DISCUSSIONS 

Our motivation to study the PN kinematics of IC 10 is based 
on the fact tha t some PN candidat es were previously found 
in this galaxy (jMagrini et al.ll2003l ), and, moreover, located 
at large galactocentric distances (around 3 kpc, see Figure[l} 
from its centre, PN 1 to 3 in Table 2. The latter can represent 
a trace of past tidal interactions and/or the connection of 
the old-intermediate stellar population with the H I envelope 
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extending 7 times farther than the optical diamet er of the 
galaxy (Huchtmcicr 1 9791 : iHuchtmeier et al.lll98ll : see also 
Meatheringham et al. 1988 for the LCM). So far, no stellar 
kinematics (of individual objects) was available (for any kind 
of stars) in ICIO, thus the link between stars and gas was 
missing, and the hypothesis that the H I envelope was a later 
acquisition, rotating in the opposite direction from that of 
the inner parts of the galaxy (WM98), needed supporting 
evidence. 

Previous to the present work there have been no kine- 
matic studies of the PN population in ICIO. Although ICIO 
is a dwarf galaxy, where not many PNe can be expected to be 
found (because of its low luminosity) , we were able to collect 
a sample of 35 PNe. A small sample in absolute terms, but 
the largest detected to date in this kind of galaxy, with the 
only exception of the Magellanic clouds (see Table 1 of Ma- 
grini 2006, with the numbe r of PNe per LG galaxy, an d also 
iMagrini fc Goncalves|[2009l and lGoncalves et al.|[20]l for an 
update). It is particularly relevant that PNe form a popu- 
lation with intermediate age between the H I and massive 
young stars and the old population of red giant stars, and 
therefore a study of their kinematics should cast light on the 
problem of the different rotation patterns of the inner disk 
and its surrounding outer gas (see, e.g., SS89). 

The presence of the PN population in the outskirts of 
ICIO implies that there was active star formation in the 
outer envelope in the past, several billions of years ago. Un- 
fortunately, without a deep spectroscopic study of the PNe 
at the outskirts of the galaxy we cannot better date the star 
formation episode that gave birth to their central stars. Such 
a study does exi st, but only for the central 5 .5x5.5 arcmin^ 
region of ICIO l|Magrini fc Goncalve j lioogl l . Following the 
latter authors, the spectroscopically observed PNe, all with 
low N/O ratio and/or with very low or absent He II, are, to 
a first approximation, "old stars". And, thus, they could be 
born in a limited period of time, 7 to 11 Gyr ago, that is, dur- 
ing the first half of the age of the Universe. This is in agree- 
ment with other indications coming from photometric stud- 
ies of th e stellar population in the outskirts of IC 10, e.g., 
those bv ISanna et al.l (|2009l ) who identified stars belonging 
to different star formation episodes (massive young stars 
(~100 Myr), old stars (~13 Gyr), and int ermediate age stars 
(~0.2-7 Gyr)). In addition. [Sanna et all JioiQ' ) found a sig- 
nificant decrease of the young stellar population when mov- 
ing from the center toward the outermost regions of ICIO. 
They detected several samples of old red giant stars (well 
fitted with isochrones of 13 Gyr) up to radial distances of 
18-23 arcmin from the galactic center, and an old star excess 
at least up to 34-42 arc min from the center. Similar result s 
were also obtained by iTikhonov fc Galazutdinorol l|2009l l. 
who detect an old stellar population coincident with the H I 
envelope. Thus we can suppose that the external PNe are 
associated with the old and intermediate-age stellar popula- 
tions, since no young population is detected in the outskirts. 
Note that we did not detect any H ll region in the FOCAS 
field A, which suggests that, presently, there is no ongoing 
star formation in the outskirts of ICIO. Moreover, being a 
starburst galaxy, with a current star formatio n ra te ranging 
firom 0.02 to 0.71 M© yr'^ (see lMated[l99i and lYin et all 
I2OIO I compilations), ICIO is und ergoing an inten se and very 



recent burst of star formation (I Yin et al. 2010l) , prob ably 



involving only its more central parts. Yin et akl (|2010l ) de- 



tailed modelling of ICIO c hemical evol ution showed, as oth- 
ers previously suspected (jMagrini fc Goncalvea .2009, ), that 
ICIO was probably formed by means of a slow gas accretion 
process with a long infall timescale of ~8 Gyr. This time 
scale is particularly significant because it corresponds to the 
epoch of the formation of the PN progenitors in the central 
regions of ICIO. If the progenitors of the PNe located in the 
outskirts of ICIO had the same age of the central ones, they 
might be the signature of an epoch of particularly active star 
formation also in the outskirts of ICIO. We might suppose 
that this important episode of star formation was driven by 
the first encounter of the H I gas cloud with the ICIO proto- 
galaxy. The conclusions by Sanna et al. (2010) are in very 
good agreement with our findings, supporting the hypoth- 
esis that the H I cloud is associated with ICIO and that it 
was hosting star formation in the past. 

Concluding, the above di scussion is in agr e ement with 
the numerical simulation by iDi Matteo et al.) (|2007l ) who 
investigated the enhancement of star formation efficiency 
in galaxy interactions and mergers using numerical simu- 
lations, and considering different types of encounters, both 
direct and retrograde. They found that retrograde encoun- 
ters produce a larger star formation efficiency than direct 
encounters. They explain this behaviour by the fact that in 
retrograde encounters tides are less efficient, allowing most 
of the initial gas to stay in the gala:xy, rather than to be 
driven outwards. This great reservoir of gas can furnish the 
fuel for an intense burst of star formation in the merging 
phase. 



6 SUMMARY 

The study described in this paper was motivated by the 
fact that a previous search for PNe in ICIO found some 
PNe located at the outskirts of this starburst galaxy (up 
to 3 kpc from the centre). Their location by itself sug- 
gests that star formation was active in the past at the very 
outer regions of the galaxy. Thus, we took advantage of 
the deep imaging and spectroscopic capabilities provided 
by FOCAS@Subaru to add a number of 29 new objects 
to th e previously kno wn (Magrini et al.|[2003l : iKniazev et al.l 
i200d : iMagrini fc Goncalves, .2009) PNe of ICIO, simultane- 
ously measuring the radial velocity of most of the galaxy's 
PNe. 

The average radial velocity of 35 of our PNe is - 
335±28 kms~^, which is, within the errors, the same as the 
recession velocity of the galaxy (-348 kms"^ NED). The 
PN radial velocities were compared with the H I gas veloc- 
ities in order to examine whether the kinematic connection 
with the underlying stellar population of ICIO could repre- 
sent the trace of past tidal interactions in the Local Group. 
We found that the H I and PN population share the same 
rotation pattern, as evidenced by the r otation curves along 
any PA between 50 and 90°, following IShostak fc Skillmanl 
( 1989, and our Fig u res 3- 5 and 8), or the velocity maps 
of IWilcots fc MiUeii (|1998| , shown in Figure 7). Moreover, 
considering that the PNe located at high galactocentric dis- 
tances in ICIO have roughly t he same age as those in the 
inner regi ons, around 8 Gyr (|Magrini fc Goncalvej l2009l : 
lYin et al.|[201Q ). the first encounter of the H I gas with the 
galaxy can be constrained to have occurred at that time. 
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Finally, our deep [O III] 5007 A FOCAS images were 
used to build up the planetary nebulae luminosity func- 
tion. From the latter we derived a ICIO distance modulus of 
24.1±0.25, which corresponds to D = 660 kpc. This dis- 
tance is in very good a greement with oth er recent ICIO 
distance dete rminations (|Sakai et al. I Il999l . 24.1±0.2; and 
iKniazev et ai.h200S . 24.30tg ig). 
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